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Abstract 

We present a study of VO2 in the rutile and monoclinic (Mi) phases by means of all-electron 
full-potential LMTO GW calculation. Full frequency dependence and the off-diagonal matrix 
elements of the self-energy are taken into account. As a result of dynamical correlation, a satellite 
structure is found above the t2g quasiparticle peak but not below, in both the rutile and monoclinic 
phases. For the monoclinic structure, the insulating state is not obtained within the usual 1-shot 
GW calculation. We perform a simplified "self-consistent" GW scheme by adding a uniform shift 
to the conduction band levels and recalculating the quasiparticle wavefunctions accordingly. An 
insulating solution with a gap of approximately 0.6 eV is obtained, in agreement with experiments. 

PACS numbers: 71.20.-b, 71.30. +h, 71.20.Bc 
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I. INTRODUCTION 



Vanadium dioxide undergoes a metal-insulator transition (MIT) at the transition tem- 
perature ~ 340 K-, and much work has been done both experimentally and theoretically to 
investigate the role of electron correlations in this transition. At high temperature VO2 is 
metallic and forms tetragonal rutile (R) structure with space group PAz/mnm, while at low 
temperature VO2 forms monoclinic (Mi) structure with space group P2i/c, and becomes 
insulating with a band gap of 0.6 eV-. The striking feature of Mi phase is the dimerization 
of vanadium atoms with the zigzag-like displacements, which leads to a doubling of the unit 
cell along the c axis. In both phases the crystal field splits vanadium t 2 g bands into a dou- 
bly degenerate band and a non-degenerate a ig band. Goodenough 3 proposed a Peierls 
(band-like) picture of the insulating state: In the metallic VO2, the a\ g band overlaps with 
the band, and both bands are partially filled. In the insulating VO2, the pairing of vana- 
dium atoms leads to the bonding- ant ibonding splitting of two a\ g bands, which separates 
the (bonding) a\ g band and e g band. The bonding a\ g band becomes filled, while the e g 
band becomes empty, and the gap opens up between them. This scenario is also supported 
by Wentzcovitch et alA-, who showed that first-principles molecular dynamics calculations 
within the local density approximation (LDA) of density functional theory can reproduce 
the structure and energy difference of these two phases in good accuracy. They failed to 
reproduce an insulating monoclinic VO2, but they attributed this failure to the problem of 
the LDA. However, since this Peierls model cannot explain the metastable insulating state 
of another monoclinic (M 2 ) structure^, where only half of vanadium atoms are dimerized, 
the importance of electron-electron correlations was emphasized^. To study the effect of 
correlations beyond the LDA, theoretical works using LDA + DMFT (dynamical mean field 
theory 9 ) have recently been donei&ii^, but a consensus has not been reached. Liebsch et 
al. performed a single-site multiband DMFT calculation with the quantum Monte Carlo 
method, but their LDA+DMFT result does not reproduce an insulating V02^. On the 
other hand, Laad et al. concluded from the LDA+DMFT calculation within the iterative 
perturbation theory (IPT), that the MIT is of the Mott-Hubbard type^i. To account for the 
spacial correlation between dimerized vanadium atoms, Biermann et al. performed cluster- 
DMFT calculations in this system and succeeded in reproducing the two phases^. In their 
result for insulating VO2 the antibonding a\ g band is strongly renormalized to form the 
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upper Hubbard band due to correlation, but the peak of the spectral function just below 
the Fermi energy is, interestingly, not an incoherent lower Hubbard band but a quasiparticle 
peak of the a\ g band, which indicates insulating VO2 has also a band-like nature. 

The above DMFT model approaches are parameter-dependent and only a few bands 
are taken into account. Hence, a first-principles description without adjustable parameters 
is highly desirable. The GW approximation (GWA)^^ has been successfully applied to 
calculations of excited state properties of a wide range of materials, thus the GWA may be 
appropriate for this problem. However, due to a system size of VO2 (monoclinic VO2 has 12 
atoms in a unit cell), a direct application of the GW method to this system has been limited. 
Continenza et al. applied a simple model GW scheme to insulating VO2 15 , and succeeded 
in reproducing the band gap. In their approach, the self-energy is approximated as a static 
non-local potential in which the experimental value of the dielectric constant is used. Very 
recently Gatti et al. performed a quasiparticle calculation with a simplified self-consistent 
GW scheme 16 -; they first carried out a self-consistent calculation within Hedin's Coulomb hole 
and screened exchange (COHSEX) approximation^ and used the resulting self-consistent 
one-particle wave functions and energies as an input for a one-shot GW calculation. Using 
this procedure they also obtained a band gap in good agreement with experiment. 

In the present work, we also perform GW calculations in order to gain better under- 
standing of electron correlation effects in VO2, with special emphasis on the influence of 
the diagonal self-energy and the role of dynamical effect in the self-energy, going beyond a 
model or static treatment of the self-energy. We find that the self-energy is in fact strongly 
energy dependent and it affects significantly the one-electron excitation spectrum. We also 
find that the off-diagonal self-energy can have a large influence on the quasiparticle band 
structure. Our result reproduces both the metallic and insulating VO2. 

II. METHOD 

In the GW method 1 - 4 ^, the self-energy E is written as the product of the one-particle 
Green's function G and the screened Coulomb interaction W as 

E(r, r', <*)) = ■£-[ G(r, r', u + u/)W(t, r', J)dJ . (1) 
2tt J 
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Here W is calculated within the random phase approximation (RPA). Quasiparticle wave- 
functions {fku} and energies {e^} satisfy the following equation 

[ - ^V 2 + ^xt(r)+^(r)]/ k ,(r) 
+ /s(r,r',^ w )/ k; ,(r')rfV = e^f^v), (2) 

where v ext and Vh are the external and Hartree potential, respectively. Usually Eq. (TSJ) is 
solved in a non-self-consistent way with two further approximations: first, the self-energy 
correction from the LDA exchange-correlation potential v^ A is assumed to be diagonal with 
respect to LDA wavefunctions, AE^(k,u) = (^ A \t(uj) - ^ c DA |-0kf A ) = ^AS w (k,w) . 
Second, the frequency dependence of the self-energy is simplified by the linearization around 
the LDA eigenenergies, 



as 



LDA^ 

LDA> 



Then the quasiparticle energies are evaluated as the first-order correction to the Kohn-Sham 
eigenvalues 

^ W = ^ A + ^AS w (k,^ A ), (4) 

where Zv v = 1 — A j s the renormalization factor. 

The reason why the above two approximations have been so successful for simple 
semiconductors^ 1 ^ is that in these systems the self-energy shows a smooth linear behavior 
around the Fermi energy, and the quasiparticle wavefunction and LDA wavefunction are 
almost identical. Since the validity of these two approximations is not clear for correlated 
system like VO2, in this work we perform calculations without using these approximations; a 
frequency dependence of the self-energy is explicitly taken into account, and the off-diagonal 
elements of the self-energy (in the Kohn-Sham basis) are included properly in the calculation. 

The quasiparticle energies are calculated by solving the following equation: 

det[(u; - e^ A )5^ -KAE^fcu)] = 0. (5) 

Here 3? means the hermitian part (i.e. = ^(A + A^)). When off-diagonal elements 

of AS(k, uj) are not negligible, we calculate the roots of Eq. (jSJ) by using the following 
linearization scheme: first we calculate AE(k, u) on uniform frequency meshes ujj = U\ + 



Au>(j — 1). Then in each region cjj < uj < ujj + \, we linearly interpolate AS F (k,w) as 

AS^(k,w) = ~ ^ AE Mi/ (k, o>j) 
u; i+ i - 

+ AS^(k,^ +1 ). (6) 

By substituting the above approximation for AS F (k, into Eq. (jSJ), the roots of Eq. (jSJ) 
are easily calculated by solving the generalized eigenvalue problem. We note that unlike the 
usual linearlization (Eq. ([3])), this scheme is exact in the limit Au — > 0. 

Our calculation is based on full-potential LMTO basis 17 . The product-basis technique is 
used--, and the frequency integral in Eq. ([1]) is numerically carried out along the imaginary 
axis with the contributions from the poles of the Green function added 21 . Details of the GW 



code are described in Ref. 



22|. Vanadium 3s, 3p electrons are treated as valence electrons, and 



151 unoccupied bands per V2O4 are used to compute G and W. To check the convergence 
of our calculation, we also performed calculations with 100 unoccupied bands per V 2 4 , and 
found that the calculated band gaps in the insulating state differ only by about 0.01 eV. For 
the sampling of the Brillouin zone, 6x6x6 (4x4x4) Monkhorst-Pack grid 23 is used for 
R (Mi) structure. In both phases, experimental lattice parameters are used 2 ^ 1 ^. 



III. RESULTS AND DISCUSSION 
A. Metallic tetragonal V0 2 

Figure Q] shows the band structure of metallic VO2 calculated within the LDA and the 
GWA. We find that the change in quasiparticle energy due to the off-diagonal self-energy 
is negligible so that we consider only the diagonal self-energy. The broad oxygen 2p band 
lies from —8 eV to —2 eV, and a separation between t 2 g bands (which extends from — 1 
eV to 2 eV) and e g bands (from 2 eV to 5 eV) is seen. This result is in agreement with 
previous calculations 26 . Compared to the LDA result, the self-energy correction makes the 
d-p separation larger by about 0.6 eV, but the bandwidth of the d bands does not show a 
significant change. We find that due to the dynamical correlation effects, at some k-points 
(not seen in Fig. [TJ the number of roots of Eq. (jSj) becomes greater than the number of LDA 
bands considered. In order to study the origin of this dynamical effect, in Fig. [2] we plot 
the frequency dependence of the diagonal self-energy A£ w (k, u) for £ 2ff bands at X and R 
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points as representative points. When the off-diagonal self-energy is neglected, the roots of 
the quasiparticle equation (Eq.(jS])) are given graphically by the intersections of two lines 
uj — £k? A an d ReAE w (k, uS). For these bands a noticeable peak is found at around 2 — 3 
eV in Re &H VV (k, u) , which produces the extra solutions of Eq.(jS]). Noting that in the GW 
approximation the imaginary part of the self-energy is related to W as^ 

unocc 

imE^(k,w > E F ) = -£ E(^k^q-k M |S qa ) 

q M a/3 

xImW^(q,o; - £ q -k M ) (^I^VM^/) 

x9(u -£ q _k M ), (7) 

occ 

ImE^(k,w < E F ) = 5]EE(^q-^l B *») 

xImW^(q, £ q _ kAt - u) (5 q/3 |^ q _ kM ^) 

x6»(£ q _ kM -w), (8) 

where {B m } is an arbitrary set of basis functions and the superscript c signifies the correla- 
tion part, this peak structure is due to the corresponding peak in IhiW^la;) or equivalently 
Ime _1 (w), where e(u) is the dielectric function. We actually find a strong peak in Ime _1 (w) 
around 2 eV, which may be interpreted as a sub-plasmon arising from strong transitions 
from the narrow a\ g band to empty states just above the Fermi level. It is interesting to 
note that the peak in ImE above the Fermi level is much stronger than the peak below the 
Fermi level. This difference in the strength of the peaks most likely originates from the 
matrix elements, as can be clearly seen in the above expression for ImE. 

In the metallic phase the orbital dependence of the self-energy for these t 2g orbitals is 
small, which reflects the weak orbital polarization in this phase. We also find that the self- 
energy shows some k-dependence; for the X point (Fig. [2](a)-(c)) and other k points that lies 
on the plane k z = 0, where k z is a reciprocal vector parallel to the c axis, the self-energy 
shows similar, or "isotropic" behavior, while for the R point (Fig. [2]^d)-(f)) and other k 
points that lie on k z = - plane, the self-energy is somewhat anisotropic and a peak around 
2 eV is more pronounced for a>i g (Fig. [2](e)) than (Fig. [2](d),(f)). 

To see the consequence of the sub-plasmon peak, we plot the diagonal spectral function 
A(k,tu) = Y,v ^l-*- m ^- £ LDA_^ s — (kZJ)\ f° r ^2g bands at X and R in Fig. [3J In calculating 
A(k, u), since the self-energy is calculated with the unshifted LDA Fermi energy, we shift 
the frequency dependence of ImE as ImE(u;) — > ImElw-AEp), where AE F = Ep W — E F DA , 



so as to reproduce the small inverse lifetime of quasiparticles around the Fermi energy. Close 
to the Fermi energy there are sharp quasiparticle peaks, and the sub-plasmon-originated peak 
at around 2 — 3 eV in the self-energy yields a weak satellite structure. The renormalization 
factor Z^ v of t 2g bands is about 0.5, which indicates strong transfer of the spectra weight 
to the incoherent part. In the DMFT calculation by Liebsch et al— and the cluster DMFT 
calculation by Biermann et a/.— the satellite structure is found also below the quasiparticle 
peak at around —1.5 eV, which is also observed in the recent photoemission spectral. That 
satellite, regarded as the lower Hubbard band, is not seen in our calculation. 

The self-energy of real systems differs from that of the Hubbard model in the following 
fashions. The energy scale of the self-energy of the Hubbard model is determined by U, 
which is of the order of a few eV, whereas the energy scale of the self-energy of the real 
systems is determined by the bare Coulomb interaction, which is typically one order of 
magnitude larger than the Hubbard U. As a result, ImS of the Hubbard model decays to 
zero after a few eV, above and below the Fermi level, whereas ImS of the real systems decays 
to zero at energies larger than the plasmon energy, which is of the order of tens of eV. 

As can be seen in Fig. |2] the presence of high-energy plasmon excitation causes ImS to 
increase to a large value outside the energy range of approximately -4 eV and 4 eV, in 
contrast to the Hubbard model in which ImS decays to zero outside this energy range. The 
decay of ImS enhances the satellite structure and may actually overemphasize the strength 



of the satellite, as discussed in Ref. 127J. Nevertheless, experimentally there appears to be 
some evidence that there exists indeed a satellite structure a few eV below the Fermi level^. 
In order to reproduce that peak, one must include the short range correlation effect beyond 
RPA; for this purpose it would be very interesting to apply methods such as GW+DMFT— 
to this system. 



B. Insulating monoclinic VO2 

Figure H] shows the band structure of insulating monoclinic VO2 calculated within the 
LDA and the GWA with only the diagonal part of the self-energy. As in the metallic phase 
we see the separation of the t2 9 and e g bands due to the octahedral crystal field and the 
downward shift of O 2p bands, but in this phase the occupied part of V 3d states is mainly 
of dig character, which results from the pairing of the vanadium atoms. To evaluate the 
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effect of the off-diagonal self-energy to the band structure, we compare the results with 
and without the off-diagonal self-energy in Fig. [51 In this calculation the matrix elements 
of the self-energy are calculated within the V d subspace. We find that inclusion of the 
O 2p bands changes the quasiparticle energies by less than 0.1 eV. In most bands the off- 
diagonal self-energy is very small and the two results are almost the same, but a noticeable 
difference is found near the Fermi energy. In the diagonal-only GW result (Fig. Efb)), an 
unusual band overlap around the A point is found, which is removed when the off-diagonal 
self-energy is included. This anomaly is attributed to the inefficiency of both the LDA 
and 1-shot GW. Around the A point the valence (ai g ) band and conduction (e£) band are 
barely separated in the LDA level (Fig. 0(a)), so these bands are "fictitiously" hybridized. 
Since the diagonal-only GW calculation does not change the character of wavef unctions, 
this hybridization cannot be removed unless the off-diagonal self-energy is included. Thus, 
the LDA wavefunctions and the quasiparticle ones obtained by including the off-diagonal 
self-energy are far from identical. 

At first sight the presence of bands crossing the Fermi level at around the A point intu- 
itively suggests that it is unlikely to open up a gap around that point. In the case of single 
band crossing the Fermi level there appears to be no choice other than an opening of a Mott 
gap. However, in the case of VO2 there are multi-bands, which still allow for an opening of a 
gap within a one-particle picture by means of rearrangement of the band occupation. This is 
indeed the case and it is noteworthy that the rearrangement of the band occupation around 
the A point is already obtained in the one-shot GW calculation, provided the off-diagonal 
self-energy is taken into account. 

We thus observe two important ingredients in gap opening in narrow band materials 
where entangled bands cross the Fermi level: First, the off-diagonal self-energy is crucial 
in " dehybridizing" the bands and rearranging the band occupation, and second, the mod- 
ification of the one-particle energies decreases the screening and hence enhances the gap 
opening 2 - 1 ^ in the course of self-consistency. This is in contrast to the conventional semi- 
conductors where the gap opening, in the case of overlapping valence and conduction bands, 
is simply affected by shifting the bands, i.e., by the second mechanism. 

We also find that due to the non-linear behavior of the self-energy, it is crucial to take 
into account the frequency dependence of the self-energy explicitly in calculating {s^f}. 
The results with and without the linearization are compared in Fig. EJ The usual linearized 
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approach (Eq. (0J), which is used in most 1-shot GW calculations, causes an error of as 
much as 0.2 eV for the a\ g band, which leads to the underestimation of the valence band 
width and hence overestimation of the band gap. Figure [7] (a-c) shows the diagonal self- 
energy for the insulating phase at V and E for three bands: bonding and antibonding a l9 and 
e* g . The self-energy for other k-points shows similar behavior. A nonlinear behavior of the 
self-energy is clearly seen from the figure, which explains the failure of usual linearization 
scheme. Owing to the lattice distortion, the self-energy for these t2 g bands shows strong 
orbital dependence, in contrast to the metallic phase; peak structures in ReS around —4 eV 
and 1 — 3 eV arising from strong transition from the narrow a± g band to empty states just 
above the Fermi level are more prominent for ai g bands. This peak yields extra solutions of 
Eq. (jSJ), which are seen as flat bands around 2.5 eV in Fig. |H The renormalization factor 
for these bands is about 0.5 — 0.6. 

In our 1-shot GW result, even if the off-diagonal self-energy is included there still exists 
a small indirect band overlap of 0.02 eV in contrast to the experimental band gap of 0.6 
eV. Furthermore, the bonding- antibonding splitting of a± g bands is less than 2 eV, notice- 
ably smaller than experiment 2 . The present calculation uses the LDA wavefunctions and 
eigenenergies to construct G and W. Thus, the failure to reproduce the experimental gap 
may be due to the poorness of the starting states, and a calculation starting from a better 
mean-field solution will be needed to test the initial state dependence. Indeed, very recently 
Gatti et al attempted such a calculation^. They first performed a static COHSEX calcu- 
lation self-consistently, and obtained an insulating phase. The result is used as the starting 
Hamiltonian of a subsequent one-shot GW calculation, which yields a final band gap in 
good agreement with experiment. They found that without modification of the LDA wave 
functions, the gap is not opened. 

Since a self-consistent GW scheme^ 2 - 1 ^ 1 ^ 1 ^ 1 ^ is not yet well-established 3 ^* 3 ^, it would be 
meaningful to study the effect of initial states 37 and self-consistency. The problem is how 
to construct a one-particle Hamiltonian whose eigenf unctions and energies well represent 
the quasiparticle wavefunctions and energies. One possible way is a so-called quasiparticle 
self-consistent scheme (QSGW) proposed in Ref. 38|, where the exchange and correlation 
part of the one-particle Hamiltonian is replaced by a static non-local exchange-correlation 
potential constructed using the following formula: 
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(9) 



with the self-energy £ calculated in the GWA. The reason for introducing such a recipe is 
that the self-energy is energy dependent. For the diagonal elements (fi = v) no ambiguities 
arise in choosing the energy but for the off-diagonal elements the choice of the energy is 
ambiguous. The formula, however, still awaits theoretical justification. 

As another choice, we propose the following procedure. The quasiparticle wavefunctions 
{^kf}) obtained as solutions of Eq. (j2j) with ImE neglected, are in general not orthogonal 
and therefore unsuitable as an input for a GW calculation. We envisage that a natural and 
physically motivated quasiparticle Hamiltonian is 



which may be viewed as a series of quasiparticle modes (QPM). This quasiparticle Hamil- 
tonian is evidently Hermitian and ambiguities in choosing the energy in the self-energy do 
not arise. Upon diagonalizing Hqp we obtain an orthonormal set of wavefunctions with 
the corresponding eigenvalues. In most cases the quasiparticle wavefunctions are almost 
orthogonal, thus the solutions of Eq. (fTDl) should be very close to the original quasiparticle 
states. 

To convince ourselves that our QPM approximation is physically sound, we compare in 
Fig. [8] the quasiparticle band structures generated from the solutions of Eq. (jSJ) and the 
QPM approximation (Eq. (flOj) ) as well as the QSGW (Eq. (jUJ)). In this calculation, the 
self-energy E is calculated with LDA eigenstates (i.e., non-self-consistent calculation). As 
can be seen the QPM band structure agree with the true quasiparticle band structure well 
within the accuracy of the calculations whereas the band structure calculated using Eq. (J9j) 
shows a noticeable deviation from the true quasiparticle band structure around the A point. 
This deviation originates from the approximated form of the off-diagonal part of the effective 
exchange- correlation potential. Thus, the quasiparticle energy can be sensitive to the choice 
of the energy in the self-energy. On the other hand, the QPM approximation is free from 
ambiguity in choosing the energy argument in the self-energy and the off-diagonal self-energy 
is properly taken into account implicitly via the quasiparticle wavefunction (Eq. (flOl) ). 

A self-consistent scheme is now at our disposal. We start by performing a GW calculation 
using the LDA band structure as an input. The resulting self-energy is then used to construct 



if QP = £i* kl ,>e 




(10) 
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the quasiparticle wavefunctions and energies by solving Eq. ((Sj). A new set of wavefunctions 
and energies are obtained within the QPM approximation in Eq. fflOl) and used to perform 
the next GW calculation. The iteration continues until the quasiparticle wavefunctions and 
energies converge. For the case of VO2 that we are considering, this is a very cumbersome 
task due to the large system size. We have observed, however, the following result: Starting 
from the metallic LDA band structure, the quasiparticle band structure obtained by solving 
Eq. (jSJ) acquires a (direct) gap already in the first iteration. The self-energy is assumed to be 
diagonal except within the V d subspace. It is important to note that the gap around the A 
point is not opened up if the off-diagonal elements of the self-energy within the tzg subspace 
are not taken into account. This suggests that the rearrangement of the wavefunctions 
already takes place in the first iteration and after the first iteration we may then keep the 
same wavefunctions. Thus in the subsequent iterations we simulate the self-consistency by 
uniformly shifting the conduction band quasiparticle energies as E^ c — > E kc + A. We tested 
this simplified self-consistent scheme for simple semiconductors and found that it yields 
a result very close to that of the full QPM self-consistent calculation. In the calculation 
of VO2, we start the calculation from the gapped state by adding an initial shift to the 
LDA conduction levels. This process is necessary to better simulate the insulating state in 
constructing the self-energy matrix. We performed the calculations with several values for 
the initial shift, and found that the final result is not sensitive to the value of the initial 
shift. 

The result of the g function of the shift is shown in Fig. [91 Two indirect gaps are 

shown, one between points B and C and the other between D and C. The GW gaps increase 
almost linearly with the shift, and a self-consistent value of 0.6 eV is found within the QPM 
approximation. This result differs considerably from the result obtained by assuming that 
the self-energy is diagonal, i.e., the wavefunctions are given by those of the LDA. Our results 
show that the GWA within the quasiparticle concept is sensitive to the treatment of the off- 
diagonal elements of the self-energy, which modify the LDA wavefunctions. Although a 
proper definition of effective quasiparticle Hamiltonian is debatable^ 1 ^ 1 ^, the obtained 0.6 
eV gap in our calculation, which happens to be in good agreement with the experimental 
value, is an indication that our QPM scheme could furnish a suitable way of constructing a 
quasiparticle Hamiltonian. 

Our result makes us reconsider the electronic state of the paramagnetic insulating VO2; 
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in Peierls description, the bonding-antibonding separation of a± g bands is determined only 
by the hybridization between two a\ g bands, and both spin up and down states are assumed 
to occupy the same (bonding a\ g ) orbitals. However, the Coulomb interaction between 
localized d electrons invalidates this interpretation, as is clear from an example of 2-site 
Hubbard model. Thus, the paramagnetic insulating VO2 may be the result of a singlet state 
formed by two localized electrons on paired vanadium atoms, which cannot be described 
in the mean-field theory This is indeed the conclusion drawn from the work of Biermann 
et al 12 , who performed cluster LDA+DMFT calculations taking the vanadium dimer as a 
unit. However, our result along with previous calculations^ 1 ^ suggest that a one-particle 
description can still go a long way in describing the electronic structure of VO2, at least 
in describing the insulating gap. Nevertheless, we think that the Peierls picture alone may 
not be sufficient for a complete description of the electronic structure of the insulating VO2. 
The missing satellites below the Fermi level, observed in photoemission experiment, is one 
such indication. 



IV. CONCLUSION 



In conclusion, we have investigated the self-energy and band structure of the metallic 
and insulating phases of VO2 within the GWA. In agreement with the work of Gatti et 
al, we have found that the GWA is able to describe the metal-insulator transition. Our 
calculations indicate that the band gap depends sensitively on the approximation used for 
the off-diagonal elements of the self-energy. Thus, a proper way of defining the off-diagonal 
elements of the self-energy is crucial. Using the QPM approximation and the scissor operator 
we have performed a self-consistent calculation to determine the band gap, which is found to 
be 0.6 eV, in agreement with calculated value by Gatti et al albeit using different methods. 
Our results within the QPM approximation suggests that the GWA can treat correlation 
effects in VO2 and the opening of the gap supports the Peierls picture of Wentzcovitch et 
al. It does not seem to be necessary to perform calculations with broken spin symmetry 
(antiferromagnetic structure) as commonly done in order to open up a gap^. 

In both the metallic and insulating phases a satellite feature exists above the Fermi energy 
but not below, in contrast to the spectra calculated within the LDA+DMFT scheme. Since 
the experimental photoemission spectrum displays satellite features above and below the 
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Fermi level, usually interpreted as the upper and lower Hubbard bands, we attribute the 
discrepancy of the GW spectra to vertex corrections, i.e., correlations beyond the RPA, 
which are apparently captured in the LDA+DMFT. 

We have also found that the common procedure of calculating the quasiparticle energy 
by linearization of the self-energy failed badly in the case of VO2. It would seem that for 
narrow band materials it is important to take into account the full energy dependence of 
the self-energy when calculating the quasiparticle energy. 
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LDA 




FIG. 1: Band structure of metalic VO2 calculated with the LDA (solid lines) and GW with only 
the diagonal self-energy (circles). 
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FIG. 2: Diagonal self-energy AE w (k, w) for three t2 g bands of metallic VO2 at the X point ((a)- 
(c)) and R point ((d)-(f)). In (b) and (e) the wavefunctions have mainly a\ g character, and in 
(a)(c)(d)(f) the wavefunctions have mainly character. Red lines: Re AE^k, uj). Blue lines: 
Im AS I/i/ (k, uj). Black solid lines: u — e^? A - The or ig m of the horizontal axes is set to the 
renormalized Fermi energy Ep^, and E]P A is shown as vertical dash-dotted lines. 
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FIG. 3: Spectral function A(k, uS) for t2 g bands of metallic VO2 at the X point (solid line) and R 
point (dashed line). 
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FIG. 4: Band structure of insulating VO2 calculated with the LDA (solid lines) and GW including 
only a diagonal part of the self-energy (circles). The symmetry labels are the same as those in 



Ref. 



26. 
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FIG. 5: Band structure of insulating VO2 near the Fermi energy calculated with (a)LDA, (b)GW 
including only a diagonal part of the self-energy, (c)GW including the off-diagonal self-energy. 
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FIG. 6: GW band structure of insulating V02- Triangles (circles) correspond to the result 
with(without) the linearization of the self-energy. Only the diagonal self-energy is included in 
the calculation. 
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FIG. 7: Diagonal self-energy AS w (k,u) of insulating VO2 for t2 g bands in the 1-shot GW calcu- 
lation at the r point (a-c) and E point (d-f). (a) and (d) : bonding a\ g bands, (b) and (e) : 
bands, (c) and (f) : anti-bonding a\ g bands. Red lines: Re AT, uu (h, ui). Blue lines: Im AT, uu (k, to). 
Black solid lines: lo — £^ A - The origin of the horizontal axes is set to the renormalized Fermi 
energy Ep W , and E]P A is shown as vertical dash-dotted lines. 
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FIG. 8: Band structure of insulating VO2 calculated with (a) GW(circles) and QSGW(squares) 
and (b) GW (circles) and QPM(triangles). 
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FIG. 9: GW indirect gap between (a)B and C and (b) D and C as a function of the initial band 
gap with the self-energy calculated by using the wavefunctions obtained within the LDA(circles) 
and QPM(triangles). 
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